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Optimizing of channel equalizer 

FIELD 

The invention relates to estimating noise power in a radio receiver 
in order to determine channel equalizer parameters. 

5 BACKGROUND 

Radio receivers employ different channel equalizers to remove 
intersymbol interference (ISI), which is caused by linear and non-linear distor- 
tions to which a signal is subjected in a radio channel. Intersymbol interference 
occurs in band-limited channels when the pulse shape used spreads to adja- 

10 cent pulse intervals. The problem is particularly serious at high transmission 
rates in data transfer applications. There are many different types of equaliz- 
ers, such as a DFE (Decision Feedback Equalizer), an ML (Maximum Likeli- 
hood) equalizer and an MLSE (Maximum Likelihood Sequence Estimation 
Equalizer), the two latter ones being based on the Viterbi algorithm. 

15 It is widely known that the information received from equalizers 

based on the Viterbi algorithm for soft decision making in decoding must be 
weighted taking noise or interference power into account in order to enable the 
performance to be optimized. The problem is then how to estimate the noise 
power in a reliable manner. 

20 Publication US 5,199,047 discloses a method which enables re- 

ception quality to be estimated in TDMA (Time Division Multiple Access) sys- 
tems. In the method, channel equalizers are adjusted by comparing a training 
sequence stored in advance in the memory with a received training sequence. 
A training sequence is transmitted in connection with each data transmission. 

25 The publication discloses a widely known receiver structure wherein impulse 
response H(O) of a channel is determined by calculating the cross-correlation 
of received training sequence X' with sequence X stored in the memory. This 
impulse response controls a Viterbi equalizer. The publication discloses a 
method which enables the reception quality to be estimated by calculating es- 

30 timate S for a received signal 
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y. is the calculated estimate for a signal (including a training se- 
quence) transmitted without interference, and 
x'. is the received sample. 

The lower estimate S is, the higher the correlation of the estimated 
training sequence with the received signal sample. Hence, the lower estimate 
S is, the higher the likelihood that the transmitted data bits can be detected by 
the channel equalizer used. 

The publication also discloses a relative estimate, i.e. quality coef- 
ficient Q, which takes the power of the received signal into account 
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wherein quadratic values of training sequence Jf,.' or individual 
sample values x'. are summed in order to determine received signal energy. 

W 15 A receiver usually, e.g. in a GSM (Global System for Mobile Com- 

' . munications) system modification called EDGE (Enhanced Data Services for 

ry GSM Evolution), comprises prefilters before the channel equalizer. Publication 

0 US 5,199,047 does not disclose how this fact can be utilized in optimizing the 

1 g i 

? S channel equalizer. 

20 BRIEF DESCRIPTION OF THE INVENTION 

An object of the invention is thus to provide a method for optimizing 
a channel equalizer by estimating noise power in two stages, and an appara- 
tus implementing the method. This is achieved by a method for carrying out 
channel equalization in a radio receiver wherein an impulse response is esti- 

25 mated, noise power is determined by estimating a covariance matrix of the 
noise contained in a received signal before prefiltering, and tap coefficients of 
prefilters and an equalizer are calculated. The method comprises determining 
the noise power after prefiltering by estimating a noise variance, and weighting 
input signals of the channel equalizer by weighting coefficients obtained by 

30 estimating the noise variance. 

The invention also relates to a radio receiver comprising means for 
estimating an impulse response, means for determining noise power of a re- 
ceived signal by estimating a covariance matrix of the noise contained in the 
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received signal before prefiltering, and means for calculating tap coefficients of 
prefilters and a channel equalizer. The receiver comprises means for deter- 
mining the noise power after prefiltering by estimating a noise variance, and 
the receiver comprises means for weighting input signals of the channel 
5 equalizer by weighting coefficients obtained from the noise variance estima- 
tion. 

Preferred embodiments of the invention are disclosed in the de- 
pendent claims. 

The invention is based on estimating the noise power, i.e. noise 
10 variance, of a received signal not only before but also after prefiltering. 

Weighting coefficients obtained from the estimation are used for weighting an 
input signal of a channel equalizer. 

The method and system of the invention provide several advan- 
tages. By weighting the input signal of the channel equalizer, the performance 
15 of channel decoding can be improved. This is particularly advantageous if, due 
to the modulation method of the system, the performance of channel decoding 
is of considerable importance, such as in a GSM modification called EDGE. In 
addition, estimating the noise again after prefiltering enables errors occurred in 
M the prefiltering to be taken into account. 

0 20 BRIEF DESCRIPTION OF THE DRAWINGS 

□ The invention is now described in closer detail in connection with 

' * the preferred embodiments and with reference to the accompanying drawings, 

in which 

. Figure 1 illustrates an example of a telecommunication system, 
25 Figure 2 is a flow diagram showing method steps for estimating a 

noise covariance twice, and potentially unbiasing an estimate, 

Figure 3 shows an impulse response of a received signal, 
Figure 4 shows a solution for calculating channel equalizer pa- 
rameters in a receiver. 

30 DESCRIPTION OF THE EMBODIMENTS 

The invention can be applied to all wireless communication system 
receivers, in network parts, such as base transceiver stations, and in different 
subscriber terminals as well. 

Figure 1 illustrates, in a simplified manner, a digital data transfer 
35 system to which the solution of the invention can be applied. The system is 
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part of a cellular radio system comprising a base transceiver station 104 hav- 
ing a radio connection 108 and 110 to subscriber terminals 100 and 102 that 
can be fixedly positioned, located in a vehicle or portable terminals to be car- 
ried around. The transceivers of the base transceiver station are connected to 
5 an antenna unit, which is used for implementing a duplex radio connection to a 
subscriber terminal. The base transceiver station is further connected to a 
base station controller 106, which conveys the subscriber terminal connections 
to other parts of the network. In a centralized manner, the base station con- 
troller controls several base transceiver stations connected thereto. 
10 The cellular radio system may also be connected to a public 

switched telephone network, in which case a transcoder converts different 
digital speech encoding modes used between the public switched telephone 
O network and the cellular radio network into compatible ones, e.g. from the 64 

l% kbit/s mode of the fixed network into another (e.g. 1 3 kbit/s) mode of the cel- 

\Q 15 lular radio network, and vice versa. 

Yi Figure 2 is a flow diagram showing method steps for estimating a 

hg noise variance in two stages, and for weighting an input signal of a channel 

' p *' equalizer by weighting coefficients obtained from the noise variance estima- 

3 

m tion. The individual method steps of the flow diagram will be explained in 

rp 20 closer detail in connection with the description of a receiver structure. The pro- 

57p cess starts from block 200. 

Q In block 202, an impulse response is calculated. Figure 3 illustrates 

^ a measured impulse response by way of example. In a typical cellular radio 

environment, the signals between a base transceiver station and a subscriber 

25 terminal propagate taking several different routes between a transmitter and a 
receiver. This multipath propagation is mainly caused by a signal being re- 
flected from surrounding surfaces. Signals propagated via different routes ar- 
rive at the receiver at different times due to a different propagation delay. This 
applies to both transmission directions. This multipath propagation of a trans- 

30 mitted signal can be monitored at the receiver by measuring the impulse re- 
sponse of a received signal, in which the signals that have different times of 
arrival are shown as peaks proportional to their signal strength. Figure 3 illus- 
trates the measured impulse response by way of example. The horizontal axis 
300 designates time and the vertical axis 302 designates the strength of the 

35 received signal. Peak points 304, 306, 308 of the curve indicate the strongest 
components of the received signal. 
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Next, in block 204, a covariance matrix of the signal is estimated, 
the diagonal thereof providing a noise variance in a vector form, according to 
Formula 7. In block 206, tap coefficients of prefilters and a channel equalizer 
are calculated using a known method. In block 208, the noise variance is esti- 
5 mated again, according to Formula 10. Finally, in block 210, the signals sup- 
plied to the channel equalizer are weighted by weighting coefficients obtained 
by the noise estimation. Arrow 212 describes the repeatability of the method 
according to the requirements of the system standard being used, e.g. time 
slot specifically. In block 214, the level of possible biasing in the estimate is 
10 assessed in order to determine parameters according to Formula 1 1 . This step 
is not necessary but will improve the performance if the tap coefficients of the 
prefilters have been determined using an equalizer algorithm which causes 
biasing to the noise energy estimate. The process ends in block 216. 

Next, each method step will be described in closer detail by means 
3 15 of a simplified receiver structure necessary for determining the channel equal- 
izer parameters, the structure being shown in Figure 4. For illustrative reasons, 
the figure only shows receiver structure parts relevant to the description of the 
invention. 

Estimation block 400 receives the sampled signal as input, and the 
20 impulse response of each branch is estimated according to the prior art by 
cross-correlating received samples with a known sequence. A method for es- 
timating impulse responses applicable to the known systems, which is applied 
e.g. to the GSM system, utilizes a known training sequence attached to a 
burst. 16 bits of the 26-bit-long training sequence are then used for estimating 
25 each impulse response tap. The structure usually also comprises a matched 
filter to reconstruct a signal distorted in the channel to the original data stream 
at a symbol error likelihood which depends on interference factors, such as 
intersymbol interference ISI. The autocorrelation taps of the estimated impulse 
response are calculated at the matched filter. The facilities described above 
30 can be implemented in many ways, e.g. by software run in a processor or by a 
hardware configuration, such as a logic built using separate components or 
ASIC (Application Specific Integrated Circuit). 

After estimating the impulse response, the noise covariance matrix 
is calculated in block 402. According to the prior art, the covariance matrix can 
35 be estimated e.g. as follows: 
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In a linear case, a sampled signal vector can be shown in the form 
(variables in bold characters being vectors or matrixes) 



y 2 ^HjX + Wj 
y 2 =H 2 x + w 2 



(3) 



wherein 

and y 2 are sample vectors of the form [yhM" + \]...y[N - i\] T , 

when n = 0,1,...,jV — 1 , wherein n is the number of samples and T 
is a transpose , 

x is the vector to be estimated, 

w 1 and w 2 are noise vectors of the form [w[«]w[« + l]...w{N- \]] T , 

H is a known observation matrix whose dimensions are 
N x (N + h, - 1) , wherein h, is the length of the impulse response 
and wherein hQ are impulse response observation values, and 
which is of the form 



H 



7?(0) h(l) ... h(h,) 0 0 ... 0 
0 h(0) h(l) ... h(h,) 0 ... 0 



0 0 



0 h(0) h(l) 



i.e. matrix H comprises an upper triangle matrix and a lower trian- 
gle matrix whose value is 0. Matrix multiplication Hx calculates the impulse 
response and information convolution. 

Thus, the covariance of the two samples and y 2 is 

*\2= E 
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wherein Eiy^ is the expected value of and of the form 
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In Formulas (5) and (6), p designates a probability density function 
and * designates a complex conjugate. 

£(y 2 )' s obtained in a similar manner. 

The covariance can be expressed in a matrix form also in the fol- 
lowing manner: 

C = E(e.e. H ^ , wherein (6) 
H designates a complex conjugate transpose of the matrix 



e - 




j y /2- H /2 x )\ 
wherein T designates a transpose of the matrix. 



(7) 



20 According to Figure 4, there may be more sample vectors than y, 

and y 2 shown in the formulas for the sake of simplicity. The elements of the 
diagonal of the covariance matrix form the signal noise variance in the vector 
form. 

The facilities described above can be implemented in many ways, 
25 e.g. by software run in a processor or by a hardware configuration, such as a 
logic built using separate components or ASIC. 

In block 404, the tap coefficients of prefilters f 1t f 2 , etc. f n and the 
channel equalizer 412 are calculated. The output signals of blocks 400 and 
402 serve as input signals of the block. The estimated impulse response val- 
30 ues and the noise covariance matrix can be used for determining the tap coef- 
ficients of the prefilters. The prefilters may be either of FIR (Finite Impulse Re- 
sponse) or IIR (Infinite Impulse Response) type but not, however, matched 
filters. IIR filters require less parameters and less memory and calculation ca- 
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pacity than FIR filters that have an equally flat stop band, but the MR filters 
cause phase distortion. As far as the application of the invention is concerned, 
it is irrelevant which filter or method of design is selected, so these will not be 
discussed in greater detail in the present description. Different methods for 
5 designing filters are widely known in the field. An output signal 416 of block 
404, which is supplied to weighting means 410, is a modified impulse re- 
sponse. 

Several channel equalizers of different type are generally known in 
the field. In practice, the most common ones include a linear equalizer, DEF 
10 (Decision Feedback Equalizer), which is non-linear, and the Viterbi algorithm, 
which is based on an ML (Maximum Likelihood) receiver. In connection with 
the Viterbi algorithm, the equalizer optimization criterion is the sequence error 
likelihood. Conventionally, the equalizer is implemented by means of a linear 
filter of the FIR type. Such an equalizer can be optimized by applying different 
15 criteria. The error likelihood depends non-linearly on the equalizer coefficients, 
so in practice, the most common optimization criterion is an MSE (Mean- 
Square Error), i.e. error power 

2 

, wherein (8) 

J . is the error power minimum, 
min 

20 1^ is a reference signal, and 

1^ is the reference signal estimate, and 
E is the expected value. 

As far as the application of the invention is concerned, 
25 it is irrelevant which equalizer or method of optimization is selected, so these 
will not be discussed in closer detail in the present description. Different meth- 
ods for optimizing equalizers are widely known in the field. 

In block 406, the signal noise variance is calculated again after 
prefiltering. According to the prior art, the noise variance can be estimated e.g. 
30 as follows: 

After prefiltering, the signal vector can be expressed in the form 
y c = H c x + w c , wherein (9) 
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min 
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y is a sample vector of the form [^[«]>-[« + l]...y[N - 1]] , when 

n = 0,l,...,N -1 , wherein n is the number of samples and T is a 
transpose, 

x is the vector to be estimated, 

w c is a noise vector of the form [w[«]i4« + l]...w[A^-l]]^, 

H is a known observation matrix whose dimensions are 

c 

N x (N + h, - 1) , wherein h c () are impulse response observation 
values and h, is the length of the impulse response, and 
which is of the form 

h c (0) h c (\) ... h^hj) 0 0 ... 0 



i°% H = 

%a# C 



o h c (V) h c {\) ... h c (hj) o 



0 



0 0 



0 h c (0) h c {\) ... h c (h l )_ 
Thus, noise energy Ncan be estimated by using the formula 
15 N= c* w£w c *!length[v/^j , wherein (10) 



c is a constant selected by the user, which is not necessary but 
which can, if necessary, be used for e.g. scaling the system dynamics, 
O length is the length of the vector, 

20 t is the transpose of the vector, 

* is a complex conjugate, and 
/ is division. 

The functionalities described above can be implemented in many 
25 ways, e.g. by software run in a processor or by a hardware configuration, such 
as a logic built using separate components or ASIC. 

If the tap coefficients of the prefilters have been determined by us- 
ing an equalizer algorithm which causes biasing to the noise energy estimate, 
such as an MMSE-DFE (Minimum Mean-Square Equalizer - Decision Feed- 
30 back Equalizer) equalizer algorithm, the estimate is unbiased in order to im- 
prove the channel encoding performance. In block 408, the weighting coeffi- 
cients for unbiasing are calculated from the noise energy estimate as follows: 
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, wherein 
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10, and 



N is the noise energy estimate and of the form shown in Formula 



is the expected value of the signal energy after prefilter- 
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This is a solution in accordance with Figure 4. 

In formula 10 for calculating noise energy N 



£w c */fe»«rt(w c ), 



N= c* w 

constant c can be determined using Formula 1 1 , already taking the 
unbiasing of the noise energy estimate into account when calculating the 
weighting coefficients. After estimating the noise energy and assessing the 
effect of potential biasing, the output signal, i.e. the modified impulse response 
416, of block 404 and a sum signal 418 formed in an adder 414 of the prefil- 
tered sample signals are multiplied by the obtained weighting coefficients us- 
ing the weighting means 410 before the actual channel equalizer block 412. 
This gives more reliable symbol error rate values for channel decoding. 

The functionalities described above can be implemented in many 
ways, e.g. by software run in a processor or by a hardware configuration, such 
as a logic built using separate components or ASIC. 

Although the invention has been described above with reference to 
the example of the accompanying drawings, it is obvious that the invention is 
not restricted thereto but can be modified in many ways within the inventive 
idea disclosed in the attached claims. 
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